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We show that the cyclic evolution of an order parameter in a many-body system that has spontaneous symme- 
try breaking in the ground state brings about a dissipative geometric phase. This phase originates from the same 
mechanism that leads to dissipation in thermo-field dynamics studied earlier by Umezawa and his collaborators 
and appears to be a manifestation of the dissipative geometric phase introduced recently by Jain and Pati. 
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I. INTRODUCTION 

Recently Jain and Pati | l|] have shown that dissipation in 
many-body systems leads to certain geometric phases. It is 
tempting to examine more closely the conditions associated 
with the appearance of those phases. It is well known that 
the low-energy response of many-body systems with sponta- 
neously broken symmetries is dominated by gapless Nambu- 
Goldstone (NG) bosons. The fact that in Jain and Pati's ap- 
proach the geometric phase is related to the static suscepti- 
bility and hence to the low-energy response of the system is 
suggestive of a possible role played by spontaneous symmetry 
breaking. 

Many-body systems with a broken symmetry of the ground 
state have certain order parameters and it is natural to ex- 
pect that the response of a time-dependent external perturba- 
tion that couples to an order parameter will cause its change 
in time. Then, following Berry's original argument |2J], one 
might expect that a cyclic evolution in the (order) parameter 
space can result in a geometric phase. Since order parame- 
ters in many-body systems are generated by symmetry break- 
ing, no additional parametric dependence of the Hamiltonian 
is needed. 

It is also well known that the Fock space of a many-body 
system is not unique but splits into folia of inequivalent rep- 
resentations, each labeled by order parameters arising from 
spontaneous symmetry breaking |3]. This implies that any 
temporal evolution of a many-body system involving change 
of an order parameter takes place within a certain class of in- 
equivalent Fock space representations. Umezawa and his col- 
laborators working in the formalism of thermo-field dynamics 
01 , 01 have already shown that temporal evolution through 
thermally inequivalent vacuua gives rise to dissipation, and 
Vitiello 1 5] has specifically pointed out that a geometric phase 
must be associated with it. 

All this taken together provides a motivation to examine 
the conditions that lead to a dissipative geometric phase in 
a many-body system with spontaneous symmetry breaking. 
We investigated here these conditions using, as an example, 
the Stoner model of itinerant electron ferromagnetism |6]. 
This model has been extensively studied in the past |7],|8], 
and various aspects of the mechanism of spontaneous sym- 
metry breaking in it have been elucidated |9],| 10]. Recently it 
was applied in the study of isospin excitations in nuclei full . 
Since, to the best of our knowledge, there exist no previous 



attempts to relate a spontaneously broken symmetry to dissi- 
pation and the geometric phase, we consider it worthwhile to 
study such a relation for a particular case of the Stoner model. 
The results obtained here appear to be generic to systems with 
broken symmetry. 

The rest of this paper is organized as follows. 

In Section II we briefly review the main features of the 
Stoner model at zero temperature following a variational ap- 
proach which gives the equation for the order parameter (mag- 
netization) from minimization of the ground state energy of 
the system. Next, in Section III we show how the emergence 
of order in the system is related to the appearance of a gapless 
collective mode, and hence to the low-energy response in the 
system. Both sections are in the nature of an abridged review, 
and the reader is referred to the original articles of Mattuck 1 9] 
and Matsumoto et al. 1 3] for details. In Section IV we study the 
closed contour integral of magnetization (the order parameter) 
as a function of an the applied external field (time dependent 
cyclic perturbation) and show that it is non-zero if the system 
is in the symmetry broken phase, and zero otherwise. That es- 
sentially gives rise to the dissipative geometric phase of Jain 
and Pati in this system and leads to the well known hysteresis 
curve in magnetic systems 1 13]. Finally, we extend our results 
to finite temperature. 



II. STONER MODEL OF FERROMAGNETISM 

The Hamiltonian of a non-relativistic many-body theory of 
interacting fermions is generic ally given by 



— ft 2 V 2 

H = I d 3 xy\x)(^—)y(x) 



i / d 3 xdW(x)^(y)V(x-y)y(yMx). (1) 



Here \|/'s are Pauli spinors and the two-body interaction po- 
tential is SU(2)-spin invariant 



V(x-y) = V Qx-y\) + Vi(\x-y\)ai-a 2 . 



(2) 



In the Stoner model this SU(2)-spin invariant electron- 
electron interaction is treated schematically as a contact-type 
interaction with a fixed strength. 



V(x~y)=X5 (3 \x-y). 



(3) 



2 



The interaction Hamiltonian reduces then to 

V = t/ 3 x\|/|(x)\|/j(x)\|/|(ji:)i|/ T (x). (4) 

The densities of spin-up and spin-down electrons are related 
to their respective Fermi momenta as 



(4) 3 



While the total electron density n remains unchanged, 



k l (4) 3 (4) 3 
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the interaction may change the individual densities of spin-up 
and spin-down electrons, giving non-zero magnetization M 



The value of A can be determined by minimizing the ground 
state energy. 

We arrive at the Stoner equation, which determines whether 
or not the system undergoes spontaneous magnetization |9] 



Xn 



— [(1 +A) 2 / 3 - (1 -A) 2 / 3 1 



(16) 



In FIG. 1, we display the plot of magnetization as a function 
of x = Xn/(2Zf). For x < 2/3, the system does not develop 



M = n\ — n\. 



(7) 



It is also useful to define the Fermi energy and momentum of 
unpolarized electron gas by 



h 2 k 2 F h 2 (3n 2 n) 2 / 3 



2m 



2m 



(8) 



We now evaluate the ground state energy of the system. The 
kinetic term is given by 



(T) 



n 2 



20m%- 



■(671 n | 



2-^5/3. 



H 2 



20m%- 



2„ , ^5/3 



(9) 



Defining a dimensionless parameter A, describing spin polar- 
ization (magnetization) by 



A = 



(10) 



one can express the nt, «| and M as functions of A 

n T = ^(l+A), Bi =^(l-A), M = nA. (11) 

Thus one arrives at 

<r) = l» eF [(l+A) 5 /3 + (i_A)n (12) 

The potential energy contribution to the ground state energy 
in Hartree-Fock appoximation is written as 

(v)hf = Mvt(*M(*)}(Yi WviW) 



(13) 



The ground state energy per electron in the Hartree-Fock ap- 
proximation (in units of Ef) is 

E(A) = 1[(1 + A) 5 ^ + (i _ A )V3] + ^(i _ A 2). (14) 



FIG. 1 : Spontaneous magnetization as a function of the density and 
coupling constant| 9] 

any spin polarization, i.e., it remains nonmagnetic. For 2/3 < 
x < 2~'/ 3 the system becomes partially magnetized, and for 
X>2-V3 k is fully magnetized. 




FIG. 2: Ground state energy of the Stoner model in Hartree-Fock 
-approximation as a function of magnetization|9] 



To summarize, in the Stoner model when a certain combi- 
nation of the density and the strength of the interaction reaches 
a critical value, the system of electrons displays a characteris- 
tic ground state instability with respect to spontaneous mag- 
netization. (In this and following figures M is used to label the 
axes, while the actual scale is given in terms of the dimension- 
less A.) 



in. MAGNETIZATION AND LOW-ENERGY SPIN 
RESPONSE 

As is well known, the appearance of an order parameter 
in a quantum many-body system implies the existence of a 
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bosonic degree of freedom in the quasiparticle Fock space of 
this system. This general statement is known as the Goldstone 
theorem ffl . 

Important as it is, Goldstone's theorem turns out to be just 
one identity out of a chain of the Ward-Takahashi| 15], 1 16] 
(WT) identities between the exact Green's functions of the 
system. These identities show how the original symmetry of 
the dynamics of the Heisenberg fields is manifested in the 
properties of the quasiparticles y|. Here we recapitulate the 
derivation of these identities according to the paper of Mat- 
sumoto et al. llOB . 

In our non-relativistic field theory of interacting fermions 
with an SI/(2)-spin invariant Lagrangian density, we add 
an external magnetic field, h, that couples to magnetization 
breaking explicitly the SI/ (2) symmetry down to 1/(1). 

L h = L(x) +h/j B ^(x)03\\f(x). (17) 
Next, we perform a global phase transformation on the 



fermionic fields 



\\t(x) — * exp[/a • a]\\t(x), 



(18) 



which we regard as a change of variables in the generating 
functional of the Green's function of our model 

W[r| t ,r|]/ 1 = —J DyDy^ exp{/ J d 4 x[L h (x) +r|fy + V t r|]}> 

(19) 

where N = W[0,0]. One then notes that, because the func- 
tional integral does not depend on the integration variables, 
the a-dependence of W/, is fictitious, and hence one should 
impose the following requirement 



dWf, 
da 



0. 



(20) 



From this condition follow the basic Ward-Takahashi identi- 
ties 



J ^x^ (x)a 3 w(x) -Y (x)a 3 y\(x))T ] t^h 
d 4 x(r\ f (x)o+y(x) - \|/ f (x)a+r\ (x) )^. n ,h 
c/ 4 x(r| + (x)c_y(x) - ^(x)o-T\(x)) n t^ h 

I — 



-2fi B h J d 4 x(a + (x)) nf ^ h 
+2fi B h / d A x{a-{x))tf nh , 



(21) 
(22) 
(23) 



where 0±(x) = (x)o±\]f(x) is the composite boson field 
and a± = G\ + i02- By taking appropriate functional deriva- 
tives of the basic Ward-Takahashi identities with respect to the 



sources, r|,r|^ and then setting these sources to zero, one ob- 
tains a chain of identities between the exact Green's functions 
in the theory. The lowest order identities in this chain are 



M h {x) = (Y(x)a 3 \\i{x)) h = 

<VT(z)v{(y))*-(vi(z)vI(y)>* = 

The first identity shows that the difference in the densities be- 
tween the spin-up and the spin-down electrons, i.e. the order 
parameter induced by the external magnetic field h, is propor- 
tional to the space integral of the propagator of the composite 
boson field. This boson field is the magnon (spin-wave). It de- 
scribes the transverse oscillations of the order parameter. The 
second Ward-Takahashi identity relates the exact propagators 
of spin-up and down electrons to their effective coupling to the 
spin-wave given by the vertex function on the r.h.s of the iden- 
tity. From our discussion in Section II we already know that 
for certain ranges of the values of the coupling constant and 
the density, non-zero magnetization may appear even in the 



-2fi B h J d A y{c + (y)C-(x)) h (24) 

-2fi B h J d*xtyi(z)y\(y)a-(x)) h . (25) 

I 

absence of an external field. This means that the l.h.s of the 
first WT-identity may be non-zero in the limit h — » 0, which 
in turn forces the r.h.s to be non-zero in that limit. In the mo- 
mentum representation this identity is expressed as 

M= — lim 2 MB h ( X( f I . ) , (26) 

where %{q 2 ) is the boson (magnon) spectral function, and 
co(q) represents its dispersion relation co = co(q). This iden- 
tity implies that in the long wavelength limit, we should have 
limq^o co(q) = and lim q ^o %(0) = M. In other words a non- 
zero magnetization in the absence of external field requires the 
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existence of a gapless, NG boson field. This in essence is the 
Goldstone theorem in the context of non-relativistic quantum 
field theory. This NG boson is physically identified with spin 
waves or magnons in this system. Its propagator has a pole at 
zero momentum which gives the dominant, low-energy spin 
response of the system. 



IV. SPONTANEOUS SYMMETRY BREAKING, 
DISSIPATION AND GEOMETRIC PHASE 

For a finite magnetic field h ^ 0, the pole contribution to 
magnetization is given in the momentum representation by 



M(h) = -2p B h 



X(q 2 ,h) 



qo-G)(q,ti) + i£ 



(27) 



<70=0,q=0 



with lim q ^ Co(q,/i) = 2p B h, and %{0,h) = M(h 



This means 

that in an external field the magnon becomes massive. This is 
a common property of pseudo NG bosons 1 3]. 

Let us assume now that the external field h is changing with 
time. Integrating both sides of equation \21\ one gets 



(28) 



M{h)dh = / %(0,h)dh. 



This essentially agrees with the result of Jain and Pati Q|, if 
one recalls that in the linear response theory the imaginary 
part of the response function (the susceptibility) corresponds 
to absorption fl2ll . In our case all the response at q = is due 
to a simple pole in the response function 



X(0,h) = -Im 



' dqo 
2k 



qo — (d(q,h)+i£ 



(29) 



q=0 



However there is also an important difference. In our ap- 
poach it suffices to couple a one-dimensional external field to 
the system, the second degree of freedom being provided by 
the order parameter itself. Thus no rotation of the external 
field is needed to generate the geometric phase. This result is 
satisfactory, since hysteresis in ferromagnets is observed for 
a magnetic source whose magnitude is varying but whose di- 
rection is fixed. 

Now we show that, when the system is in a broken symme- 
try ground state and the integral (I28t is evaluated for a certain 
closed contour in h —M plane, that integral is non- vanishing. 
In order to find the functional relation of magnetization to the 
external field M{h), we need to evaluate the Green's functions 
of quasiparticles in the system. Here we calculate the Green's 
functions approximately by truncating the (infinite) chain of 
WT-identities and solving the resulting system of equations. 



The lowest order approximation is achieved by truncating the 
WT-chain of identities at the level of the vertex-propagator 
identity 



= 2p B h-MT(p,-p,0). (30) 



If one approximates the vertex with a constant A, one can 
determine the approximate boson and fermion propagators 
consistent with the WT-identities. This corresponds to the 
Hartree-Fock (HF) approximation for the fermions and the 
random-phase-approximation (RPA) for the boson (magnon) 
field. The approximate fermionic propagators are 



Sf(p) 

s" f (p) 



l 



po - (Ep - n) + 1 (2p B h - AM) + ie 



po - (e p - p.) - \ (2/j B h - AM) + z'e ' 



(31) 



(32) 



where p is the chemical potenitial of the interacting system, 
and £ p = p 2 /(2m). Note that this structure of the spin-up 
and down electron propagators follows from the spontaneous 
St/ (2) — > f/(l) symmetry breaking. This structure is observed 
in experiments as a characteristic splitting in the dispersion re- 
lations of the spin-up and down electrons 11711 . 

Having found the Green's functions of the single-particle 
states in HF-approximatation, we now can estimate the densi- 
ties of spin-up and down electrons. For the majority spin rtf 
we get 

&(-e p +n - i (2/j B h - AM) ) ^ 
-^ m ^- l -(2^ B h-'kM)) V2 - 



(33) 



and similarly for the minority spin n\\ 



(34) 



One needs to determine p., and M, subject to the following 
constraints: 



n = n-f -j-ni 
M = n-f — «j. 



(35) 



Eqs. l|33}-!|35) lead to 



J 



M 



ffji-i (2p B h - AM)) 3 /2 + + l -{2p B h - AM)) 3 /2 
f^{p- 1 - (2p B h AM)) 3 / 2 - (p + l -{2p B h - AM)) 3 / 2 



(36) 
(37) 
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Now using the definitions in eq.l^J these conditions can be 
rewritten as 



1+A = ^(p-^Bh-XM)) 3 ? 2 



(38) 



1-A = -^^+ l -(2^Bh-lM)fl 2 . (39) 

We thus arrive at an equation linking the magnetization to the 
applied field 



(l-A) 2 / 3 -(l+A) 2 / 3 = 2 ^' XwA 



(40) 



Note that for h = it reduces to the previously obtained vari- 
ational condition that determines spontaneous magnetization. 
Thus ea. (l40> represents a generalization of ea. (!16> for a case 
of non zero external field h. One can also obtain this condition 
directly by fixing h and minimizing the ground state energy in 
HF-approximation with respect to A 



(A, A) 
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e F «((l+A) 5 / 3 + (l-A)V3) 
Xn 2 

+— (l-A 2 )+/i B n/zA. 



(41) 



The dependence of magnetization on h is illustrated in FIG. 3 

h 



FIG. 3: Driven magnetization M(h) in a symmetric ground state 
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FIG. 4: Driven magnetization M(h) in an asymmetric ground state 

and FIG.4 for the two cases: a) the system does not display 
spontaneous magnetization in its ground state and b) the sys- 
tem gets spontaneously magnetized in the absence of external 
field i.e. when 2/3 < x < 2" 1 / 3 . 




FIG. 5: Ground state energy along the valley of stability as a function 
of magnetizarion 



One can relate the characteristic s-shape of the curve M(h) 
to the form of the meandering symmetric valley of the min- 
imum ground state energy. Although the bed of this valley 
remains a local minimum for fixed h, this bed does not stay at 
a constant level but rises and falls as the applied field changes. 
One can express h through M from eq.j40>. and substitute that 
expression in E(h(M),M). The resulting level of the valley 
as a function of M is given in FIG. 5. We see that at the 
points on the hysteresis curve where M jumps (a — > b and 
c — > d) the destination points have lower energy than the start- 
ing points. This indicates that directional tunneling (in reality 
also thermally assisted) between these states may start ahead 
of the turning points, as often observed in the experimental 
hysteretic curves. 

One can see from these pictures that a system that displays 
an instability of the ground state with respect to spontaneous 
magnetization, behaves very differently from a system that 
does not have such instability. The former displays a hys- 
terisis loop when the external field is varied in a cycle, while 
the latter does not. The area of the hysterisis loop gives the 
dissipated magnetic energy. In other words we have proven 
that 



M(h)dh = & %(0,h)dh^0. 



(42) 



This area gives the cyclic Jain-Pati dissipative geometric 
phase in our model. The result can easily be extended to fi- 
nite temparature. The dominant part of the low-energy spin 
response at finite temperature again comes from the magnon 
pole at 2unh . bu t now its spectral function is temperature de- 
pendent 1 3 a , on , o 



M(h,T)=Im — 



' dqp 
2k 



l{q 2 XT) 



(43) 



q=0 



This suggests that M(h,T) can be determined from a finite 
temperature extention of ea. (l40> using the standard technique 
described in 11911 . We first define 



t 1 

el = e p + -(2 f i B h-'kM) 



p --{2fi B h-XM). 



(44) 
(45) 
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Then the spin-up electron density is 

„ T = -i& a ff*IL( S " F(Po > P) )e-^ 
n-^oj J {2%f ^l+expp^-ej,) J 

(2m) 2 / 3 r~ 

4ti 2 Jo l + exp(P0u-e-i(2^-AM)) 

(2m) 2 /\,^ 3 /2 



dz 



4k 2 



-(kT) 



2 1 



1 +exp(a(M,/;) -; 



c/jc, (46) 



with 



P = l/W x = Pe a(M,A) = P0i-|is/iH->Af/2). 
Similarly the spin-down electron density is 



(2m) 2 / 3 



A/ 2 



o 1 +exp(y(M,/i) -x 



■dbe, (47) 



with 



Defining 



y(M,h) = ^+fiBh-XM/2). 



poo 

F(U) = / r 
Jo 1 



A/ 2 



- exp(a(M, h) - x 



-dx, 



(48) 



(49) 



one finds 



M=^^(kT)y 2 (F(a(M,h))-F(y(M,h))) 

« = ^^(kT) 3 / 2 (F(a(M,h))+F(y(M,h))), (50) 



from which it follows 

'(F(a(M,h))-F(y(M,h)) 



M = n 



(F(a(M,h))+F(j(M,h)) 



(51) 



One can in principle determine both M(h,T) and p(h,T) from 
these equations by solving them numerically. Then to com- 
pute the finite-temperature Berry's phase one simply needs to 
integrate the finite-temperature order parameter along the hys- 
teresis loop. There is nothing surprising about this result. We 



have simply come to a reinterpretation of the first law of ther- 
modynamics as a balance of dynamic and geometric phases in 
a many-body system 



TdS = dU + pdV — HdB 



(52) 



As often seen before, the geometric phase gives us a new per- 
spective on a classical result. 



V. CONCLUSIONS 

We have shown that when a quantum many-body system 
with a spontaneous symmetry breaking in the ground state 
evolves in time in such a way that an order parameter com- 
pletes a cycle, one can associate a certain geometric phase 
with the energy dissipated in this cycle. We have examined 
the Stoner model of itinerant ferromagnetism and have shown 
that when the system has ferromagnetic instability, i.e., spon- 
taneous magnetization occurs in the absence of applied ex- 
ternal field, then hysteresis arises with respect to a changing 
external field and the magnetization process is irreversible. 
Thus, spontaneous symmetry breaking seems to have bear- 
ing upon another important concept in physics - reversibility. 
The measure of irreversibility is given by the dissipation, i.e., 
the dissipative geometric phase. Our preliminary investiga- 
tion has confirmed that a similar mechanism is responsible for 
the hysteresis observed in other many-body systems display- 
ing spontaneously broken symmetry in the ground state. We 
plan to address this issue in a forthcoming publication. 
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